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INTROIXJCTION 


CONGRATS  (CCNtinuous  Gradient  Ray  Tracing  System)  is  an  integrated 
collection  of  ray  tracing  programs  assigned  to  construct  ray  diagrams, 
to  'determine  multipath  eigenrays,  to  model  acoustic  propagation,  to 
compute  boundary  and  volume  reverberation,  and  to  solve  numerous 
additional  problems  in  underwater  acoustics.'  The  CONGRATS  series  has 
been  documented  in  references  (a),  (b),  and  (c). 


Bottom  reflection  loss  is  often  an  important  factor  in  acoustic 
field  predictions.  The  foundational  CONGRATS  programs,  S0990  and 
S0991,  require  am  input  table  of  bottom  loss  in  decibels  as  a function 
of  grazing  angle  in  degrees.  This  data  is  often  difficult  to  obtain 
directly.  Hence,  various  expressions  have  been  developed  to  model 
bottom  loss  versus  grazing  angle.  v 

If  a predetermined  table  is  not  supplied,  bottom  loss  values  may 
be  computed  automatically  by  one  of  four  models  which  have  been  added 
to  the  CONGRATS  programs,  5099^  and.  £0991*  The  first  model,  which 
depends  on  source  frequency,  is  a least- squares  fit  of  a fourth  degree 
polynomial  surface  to  modified  AKOS  bottom  loss  data.  The  model  was 
developed  by  Dr.  E.  S.  Eby  (see  reference  (d)).  The  second  model, 
supplied  to  the  authors  by  E.  M.  Podeszva,  is  an  approximation  to  aver- 
age bottom  losses  in  KGS  Acoustic  Provinces  1,  2,  3 , *+,  and  ?.  The  third 
and  fourth  models  are  empirical  bottom  reflection  loss  expressions, 
depending  on  both  source  frequency  and  bottom  porosity,  developed  by 
N2L  (described  in  refe^dnce  (e))and  NUNC  (see  reference  (f)),  respec- 
tively. This  memorandum  will  describe  the  four  models,  make  simple 
comparisons,  describe  the  input  options  needed  to  use  the  various 
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models,  and  provide  a listing  of  the  revised  Subroutines  INRJT  and 
BLOSS  and  the  new  Subroutine  BTMLOS. 

DISCUSSION  OF  MODELS 


BOTTOM  LOSS  TABLES 

When  an  acoustic  ray  reflects  from  the  ocean  bottom,  the  intensity 
of  the  ray  diminishes.  Tie  amount  of  power  lost  depends  on  the  frequency 
of  the  ray,  the  angle  with  which  the  ray  intersects  the  bottom,  and  the 
characteristics  of  the  bottom.  If  the  bottom  characteristics  are 
assumed  to  be  constant  over  the  area  of  interest  and  if  the  frequency 
of  each  ray  is  assumed  to  be  the  source  frequency,  the  bottom  loss 
depends  on  grazing  angle  alone.  In  order  to  calculate  bottom  less  for 
each  ray,  the  CONGRATS  programs  require  a table  of  bottom  loss  values 
for  grazing  angles  varying  from  0°  to  Q0C  in  increments  of  1°.  Linear 
interpolations  are  performed  between  these  values  in  order  to  compute 
bottom  loss  for  other  grazing  angles.  This  table  may  be  supplied  by 
the  user  on  eleven  computer  cards.  The  first  card  contains  the  words 
"BOTTOM"  and  "LOSS"  starting  in  columns  1 and  11,  respectively  (see 
Fig".  1(a)).  Ten  cards  containing  the  91  necessary  values  of  bottom 
loss  in  decibels  must  follow  the  BOTTOM  LOSS  card  (see  reference  (a)). 

The  use  of  this  table  allows  the  bottom  loss  curve  to  be  completely 
arbitrary. 

AMOS  BOTTOM  LOSS 

If  the  user  is  not  familiar  with  the  bottom  characteristics  but 
does  know  the  frequency  at  which  he  is  operating,  the  bottom  loss  table 
can  be  generated  by  the  least- squares  fit  to  modified  AMOS  data  described 
in  reference  (d).  This  fit  is  good  only  for  grazing  angles  less  than 
60°,  however.  Hence,  an  approximation  to  bottom  loss  data  (presented 
in  reference  (d))  was  developed  by  the  authors  for  grazing  angles 
steeper  than  60°.  The  combination  of  this  high-angle  approximation 
and  the  least-squares  fit  was  programmed  for  the  CONGRATS  series.  The 
equations  used  are  shown  in  Thble  1.  In  order  to  implement  the  use  of 
this  model,  the  card  snown  in  Fig.  1(b)  must  be  used  instead  of  the 
BOTTOM  LOSS  set  described  in  reference  (a).  It  should  be  noted  that  a 
FREQUENCY  card  (see  reference  (a)),  such  as  the  one  in  Fig.  2,  must 
precede  this  BOTTOM  LOSS  card.  Figure  3 is  a graph  of  the  approximated 
AMOS  bottom  loss  versus  grazing  angle  for  frequencies  of  250,  500,  1000, 
2000,  and  kOOO  Hz.  The  general  trend  shown  in  the  plot  is  that  bottom 
loss  increases  with  frequency. 
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MGS  BOTTOM  LOSS 

If  the  bottom  losses  in  the  operational  area  of  interest  are 
knovn  to  be  similar  to  those  of  one  of  the  MGS  Acoustic  Provinces,  an 
approximation  to  the  average  bottom  losses  in  Province  1,  2,  3 , 4 , or  5 
can  be  supplied  by  the  CONGRATS  programs.  The  mathematical  expressions 
used  as  the  approximations  are  shown  in  Table  2.  These  formulas  depend 
only  on  grazing  angle,  so  bottom  loss  tables  can  be  supplied  by  the  pro- 
grams with  no  additional  data.  The  approximate  bottom  losses  of  MGS 
Acoustic  Province  1,  2,  3#  4,  or  5 will  be  computed  if  the  BOTTOM  LOSS 
card  shown  in  Fig.  1(c),  (d),  (e),  (f),  or  (g),  respectively,  appears 
in  the  data  deck.  The  bottom  losses  of  the  five  MGS  Acoustic  Provinces 
are  plotted  in  Fig.  4.  For  typical  sonars  in  use  today.  Provinces  1 
and  2 provide  excellent  to  good  bottom  reflectivity.  Province  3 provides 
a fair  bottom  bounce  area,  and  Provinces  4 and  5 are  marginal  to 
unusable. 

NEL  BOTTOM  LOSS 

The  composition  of  the  bottom,  which  can  be  usefully  classified  by 
the,  single  parameter  porosity,  is  an  important  factor  in  determining 
bottom  reflectivity.  Table  3 (reproduced  from  reference  (f))  presents 
porosity  as  a function  of  bottom  composition,  with  the  bottom  varying 
from  fine  clay  (high  bottom  loss)  to  very  coarse  sand  (high  bottom 
reflectivity).  The  bottom  loss  model  developed  by  the  Naval  Electronics 
Laboratory  (and  described  in  reference  (e))  depends  on  both  source 
frequency  end  bottom  porosity,  as  well  as  grazing  angle.  The  equation 
used  is  presented  in  Table  4.  Figure  5 is  a plot  of  EEL  bottom  loss 
versus  grazing  angle  for  frequencies  of  250,  500,  1000,  2000,  and 
4000  Hz  and  a constant  porosity  of  0.60.  The  bottom  loss  increases  with 
frequency  in  a fashion  similar  to  that  shown  in  the  AMGS  approximation. 
Figure  6 is  a plot  of  MEL  bott  an  less  versus  grazing  angle  for 
porosities  varying  from  0.30  to  0.90  in  steps  of  0.15  and  a constant 
frequency  of  1 kHz.  As  bottom  porosity  increases,  bottom  loss  also 
increases.  Figure  1(h)  presents  the  data  card  required  to  cause  the 
NEL  bottom  loss  equation  to  be  used  by  the  CONGRATS  series  for  a 
porosity  of  0.60  and  a previously  determined  frequency. 

NUWC  BOTTOM  LOSS 

The  Naval  Undersea  Warfare  Center  has  developed  a newer  bottom 
loss  expression,  based  primarily  on  the  AMOS  results,  which  also  depends 
on  both  source  frequency  and  bottom  porosity  (see  reference  (f)).  The 
equation  is  presented  in  Table  4,  below  the  equation  used  by  MEL.  The 
frequency  factor,  ABTML0S,  is  a function  consisting  of  straight  line 
segments  joining  eight  specified  points.  The  eight  points  are  tab- 
ulated in  Fig.  7 and  a plot  of  the  function  is  drawn.  The  NUWC  bottom 
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loss  increases  with  frequency  2nd  porosity  (shown  in  Figs.  8 and  9, 
respectively)  as  did  the  NEL  bottom  loss.  If  uhe  card  shown  in  Fig  l(i) 
appears  in  the  data  deck,  the  NUWC  bottom  loss  model  is  used  by  the 
CONGRATS  programs  for  a porosity  of  0.40  and  a previously  inputted 
frequency. 

COMPARISON 

For  a frequency  of  3*5  kHz  the  bottom  loss  curves  derived  from 
the  AMCS,  NEL,  and  NUWC  models  (with  a porosity  of  0.40  used  in  the 
NEL  and  NUWC  approximations)  are  plotted  in  Fig.  10.  If  source  fre- 
quency is  the  only  information  available,  the  AMCS  model  produces  a 
reasonable  approximation  to  a characteristic  bottom  loss  curve.  If  a 
small,  (large)  bottom  porosity  (and  hence  high  (low)  bottom  reflectivity) 
is  used  in  the  NEL  or  NUWC  equations,  the  bottom  loss  curve  produced, 
at  any  given  frequency,  is  lover  (higher)  than  the  curve  produced  by 
the  “MGS  approximation,  which  was  developed  for  an  average  bottom. 

This  phenomenon  is  demonstrated  for  a low  porosity  in  Fig.  10.  For  a 
given  frequency  and  porosity,  the  NEL  model  usually  gives  higher 
bottom  loss  curves  than  the  NUWC  model,  especially  for  smaller  grazing 
angles.  The  general  shape  and  relationship  of  the  two  curve s is  shown 
in  Fig.  10.  If  no  other  criterion  is  available  for  choosing  between 
the  NEL  and  NUWC  bottom  loss  expressions,  it  should  be  noted  that  the 
NUWC  model  was  developed  more  recently  and  is  more  widely  accepted. 

If  the  operating  area  has  been  previously  classified  as  an  liGS  Acoustic 
Province,  the  CONGRATS  approximations  to  the  average  bottom  losses  in 
Provinces  1-5  can  be  considered  quite  reliable. 

DISCUSSION  OF  INPUT 

The  choice  of  bottom  loss  model  depends  on  the  units  field  (the 
ten-column  alphanumeric  field  starting  in  column  21  — see  reference  (a)) 
of  the  BOTTOM  LOSS  card  appearing  in  the  CONGRATS  data  deck.  A list 
of  the  nine  words  that  may  be  used  in  the  units  field,  a brief  des- 
cription of  the  model  implemented  by  each,  and  a inference  to  a sample 
card  are  shown  in  Table  5.  If  the  AMCS,  NEL,  or  NUWC  model  is  used, 
a FREQUENCY  card  (see  reference  (a)),  such  as  the  one  shown  in  Fig.  2, 
must  precede  the  BOTTOM  LOSS  card.  The  ten-column  numeric  field 
starting  in  column  31  of  the  BOTTOM  LOSS  card  contains  the  porosity 
of  the  bottom  if  the  NEL  or  NUWC  expressions  are  used. 

An  option  has  been  provided  to  input  a minimum  and/or  maximum 
bottom  loss  value  in  decibels  to  be  accepted  by  the  CONGRATS  programs. 

A bottom  loss  value  calculated  to  be  smaller  (larger)  than  the  minimum 
(maximum)  bottom  loss  is  converted  to  the  minimum  (maximum)  value.  The 
ten- column  numeric  fields  starting  in  columns  Ul  and  51  of  the  BOTTOM 
LOSS  card  contain  the  minimum  and  maximum,  respectively.  If  these 
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fields  are  left  blank  a minimum  of  0 dB  and  a maximum  of  500  d3  will  be 
assumed.  This  option  cannot  be  used  if  a bottom  loss  table  is  inputted. 


SUMMARY 


Hie  fundamental  CONGRATS  programs,  SO 990  and  S0991,  have  been 
revised  to  calculate  bottom  loss  as  a function  of  gracing  angle.  Hie 
bottom  loss  curves  have  been  computed  with  frequency,  MGS  Acoustic 
Province,  and  both  frequency  and  porosity  as  basic  parameters.  The 
bottom  loss  models  have  been  described  in  detail,  examples  compari- 
sons have  been  provided,  and  the  new  input  options  have  been  explained. 
Subroutines  INPUT  and  BLOSS  haw  been  revised  in  order  to  accommodate 
the  new  bottom  loss  models,  and  a new  subroutine.  Subroutine  BTMLOS, 
has  been  written  to  compute  the  bottom  loss  curves.  These  three 
routines  are  listed  in  the  appendix.  The  ray  plotting,  eigenray 
generation,  and  eigenray  processing  functions  of  the  CONGRATS  series 
have  not  been  altered  in  any  way  other  than  as  described  in  this 
memorandum.  It  Is  hoped  that  the  additional  bottom  loss  options  will 
be  a convenience  for  users  of  the  CONGRATS  series. 


<J  ilr  * REi  S.  COHEN 
tethematician 


TE2LDA  A.  GARRETT 
Math  Aid 
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TABLE  1 

AMOS  Equations 

Let  f ■ frequency (kHz) , 9 * Grazing  Angle (deg),  Ng  * Bottom  Loss(dB) 
If  f < 1 kHz. 


for  9 - 0°:  NB  * 0.0 


for  1°<  9 <14°:  Nfl  - 2.96026  + 1.26990F  + 1.604300  + .25257F2 

+ . 66161F0  + . O6O7202  + .02000F3  + .O9833F20 

- .O1888F02  - .0261903  - .00303F4  + .OO77SF30 

- .OO833F202  - .O125OF03  - .O131O04 
where  0 ■ (9  - 10) /5  and  F * log2  f 


for  15°  < 9 1 39°:  Ng  = 4.35024  + 1.17091F  + .468740  + .08272F2 

+ . 17115F0  - .0754702  - .00386F3  + .O1295F20 

- .O1678F02  - .0055503  - .00003F4  - .OO193F30 

- .OOO74F202  - .OO224F03  + O-04 
where  0 = (9  - 30) /5  and  F = 2 log2  f + 3 


for  40°  < 9 < 60°:  Ng  «=  4.84820  + 1.46858F  - .114570  + .11036F2 

- .O1238F0  - .OO95302  - .01000F3  - .OO375F20 

- .OO929F02  + .0083303  + 0*F4  - .OOO42F30 

. - .OO268F202  + .OO83OF03  - .0020894 

where  0 = (9  - 50)/ 5 and  F * 2 log2  f + 3 


for  61°  < © < 90°:  Ng  - 4.61430  + 1.47306F  - .0475160  + .09214F2 

- .01084F3  - .O12717F0 
where  0 = 9 - 60  and  F = 2 log,  f + 3 


NOTE;  For  61°  £ 9 £90°  Ng  can  also  be  written: 
Ng  - NB(9  - 1)  - .0S5667  - .025434F 
where  F ■ log2  f 
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TABLE  1 (Continued) 

AMOS  Equations 


If  f > 1 kHz. 


for  9 - 0°:  NB  - 0.0 


for  \°  < 9 < 14°:  NB  - 2. 96026  + 1.26990F  + 1.604300  + .25257F2 

+ . 66161F0  + .O6O7202  + .02000F3  + .O9833f20 

- .O1888F02  . .0261903  - .00303F4  + .OO778F30 

- .OO833F202  - .O125OF03  - .0531004 
where  0 = (9  - 10) /5  and  F = log2  f 


for  15°  < 9 < 39°:  NB  =*  12.25083  + 3.88808F  + 1.148020  + .16312F2 

+ .C4386F0  - . 1883502  - .02857F3  - .O82l4F20 
- .O5952F02  + .OO66703  + .01071F4  - .OO179F30 
+ .OO476F202  + .O1389F03  - .OO34804 
where  0 = (0  - 30) /5  and  F = log£  f - 1 


for  40°  < 9 < 60°:  NB  = 13.59067  + 3.74190F  - .259280  + .13266F2 

- .O3764F0  - .O881602  - .01333F3  + .O1214F20 

- . OO959F02  + .OO5OO03  - .02000F4  - .OO167F30 
+ .OO867F202  - .OO167F03  + .OOO8304 

where  0 = (9  - 50) /5  and  F = log2  f - 1 


for  61°  < 0 < 90°:  NB  - 12.77275  + 3.61490F  - .1263350  + .19162F2 
' - . OOO333F0  - .01667F3  + .O34OO20/(F  + 1) 

- .02000F4 

where  0 = 9 - 60  and  F = log2  f - 1 


NOTE:  For  61°  <0  £ 90°  NB  can  also  be  written: 

Nb  * Nb(9~-  1)  - .126002  - .000333F  + .034002/F 
where  F =•  log2  f 
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TABLE  2 
MGS  Equations 

Let  9 - Grazing  Angle (deg)  and  NB  - Bottom  Loss(dB) 


for  Province  1: 


«B 


2.2435  loge(. 12609  + 1.496) 


for  Province  2: 


Nb  = 3.4315  loge(. 10569  + 2.842) 


for  Province  3: 


nb 


= 2.4910  loge(. 88649  + 10.526) 


for  Province  4: 


NB  = 2.8377  loge(l. 87549  + 15.685) 


for  Province  5: 


NB  = 2.4036  loge(20. 57609  + 82.440) 


.vj  ***  •• ; 
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TABLE  4 

NEL  and  NUWC  Equations 


NEL  Equation: 


Nr 


6 + 22(P  - .27)  -*•  10  logl0f  J • 

TANH  ((P/ . 24) • (9/ 57. 296))  + ((1  - P/ . 24) / 1 2 . 5) • (g/90)2  J 


NUWC  Equation: 

Nb  = [3.7  + 17.5(P  - * 27) j * j^ABTMLOS ( f)  j ♦ 

f TANH ((P/. 24)- (9/57.296))  1>5/F  + ((1  - P/ . 27) /12. 5) • (e/90)2 J 


where:  Ng  = Bottom  Loss(dB) 
f = Frequency(kHz) 

P = Porosity 
9 * Grazing  Angle (deg) 


NOTE:  ABTMLOS  is  a function  of  frequency  defined 
in  Figure  7. 
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Units  Word 

Sample  Card 

Model 

Fig.  1(a) 

Inputted  Table  of  Values. 

AMOS 

Fig.  1(b) 

Least-Squares  Fit  to 

Modified  AklOS  data. 

UGS1 

Fig.  l(o) 

Average  Bottom  Losses  in 

UGS  Aooustio  Province  1. 

UGS2 

Fig.  1(d) 

Average  Bottom  Losses  in 

L5GS  Acoustic  Province  2. 

UGS3 

Fig.  1(e) 

Average  Bottom  Losses  in 

MGS  Aooustio  Province  3. 

MGS4 

Fig.  1(f) 

Average  Bottom  Losses  in' 

MGS  Acoustic  Province  4. 

MGS  5 

Fig.  1(g) 

Average  Bottom  Losses  in 

MGS  Aooustio  Province  5. 

NEL 

Fig.  1(h) 

Empirical  Expression  (NEL). 

NUWC 

-H 

H 

• 

to 

£ 

Empirical  Expression  (NU7IC). 
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Fig.  3 - AMOS  Bottom  Loco  Curves 


Pro vinos 


MGS  Bottom  Lobs  Curves 


Frequenoy 


NEL  Bottom  Loss  Curves 
(Porosity  o 0.60) 


NEL  Bottom  Loss  Curves 
(Frequency  - 1 kHz) 


TABULATED  FUNCTION 


Fig.  7 - Function  ABTOLOS 


NUWC  Bottom  Lobs  Curves 
(Porosity  ■ 0.60) 
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CHIC  weinbers 


DIMENSION  TARGET( 1010 . 3) » ANGLES ( 1010) ,S0NAR(5»2) »SURFAC(5) » 

1 SP (210.2) .BOTTOM (5) »BP(210,2) »T0lERA(5) , VP (210.2) .TP (210 .3) 

DIMENSION  GROUP ( 11 )/66HC0MMENTARGETS0NaR  SURFACBOTTOMVELOCITHERMA 
1PR0CESEND  FREQUEMAXIMU/ 

DIMENSION  COORQ ( 8 ) /48HRANGE  DEPTH  ANGLE  LOSS  TOLERAAXES  PROFILP 
1HASE  / 

DIMENSION  UNITS(2.32)/ 


1 

1.09361111E-5, 

6HCM  . 

1.09361U1E-5, 

6HCM/S 

9 

2 

2.77777777E-5, 

6HIN  . 

2.77777777E-5. 

6H1N/S 

9 

3 

0.33333333E-3, 

6HFT  » 

0.33333333E-3, 

6HFT/S 

9 

4 

1.00000000E-3, 

6HYD  . 

1.00000000E-3, 

6HYD/S 

9 

5 

1.09361111E-3, 

6HM  . 

1.0936H11E-3. 

6HM/S 

9 

6 

2.00000000E-3. 

6HF  » 

2.00000000E-3, 

6HF/S 

9 

7 

0.33333333E-0, 

6HKFT  » 

0.33333333E-0, 

6HKFT/S 

9 

8 

1.00000000E-0, 

6HKYD  . 

l.OOOUOoOOE-O, 

6HKYD/S 

9 

9 

1 . 09361 lllE-0. 

6HKM  . 

l.093bllllE-0. 

6HKM/S 

9 

T 

1.76000000E-0. 

6HMI  . 

1.76000000E-0, 

6HMI/S 

9 

1 

2.02680000E-0, 

6HN  MI  . 

2.026faOo0OE-0. 

6HKN0TS 

9 

2 

1.74532925E-2, 

6H0EG  . 

I.74532925E-2, 

6HDEG/S 

9 

3 

1.00000000E-0, 

6HRAD  . 

1.00000000E-0. 

6HRA0/S 

9 

4 

6.28318531E-0. 

6HCPS  . 

6.28316531E+3, 

6HKCPS 

9 

5 

1.00000000E-3, 

6HMS  . 

1.00000000E-0, 

6HSEC 

9 

6 

1.00000000E-0. 

6HC  . 

5.55555555E-1. 

6HFAHR 

/ 

DIMENSION  TEST (12) »FMT(12) ,DATA(2048) ,REVMAX(4) 

DIMENSION  CMMNT(12)/72H  THERE  SHOULD  BE  AT  LEAST  ONE  COMMENT  SET 
1.  / 
DATA  PRINT/6HPRINT  /.LINMAX/60/. ICMMNT/0/ 

COMMON  TARGET . ANGLES , SONAR , SU«F AC , SP  » BOTTOM , BP . TOLERA , VP . 

1 REVMAX.PROCES 


SET  INITIAL  CONDITIONS. 


ANGLES(IOIO)  = 0.0 
TARGET (1010.1)  = 0.0 
TARGET(1010.2)  = 0.0 
IhEAD  = 0 
LINES  = LlNMAX 
NPRINT  = FLO ( 33.3. PROCES) 
NTAPE1  = FLD(30. 3. PROCES) 
NTAPE2  = FLO (27 .3. PROCES) 


CHECK  THE  GROUP  CODE. 

100  READ  101.  (TEST(J) . J=1.9) 

101  FORMAT(  3(A6.A4).  3F10.5  ) 

DO  110  1=1.11 

IF ( TEST ( 1 ) . GT , GROUP ( I ) .OR.  TEST(l) .LT.GROUP(I)  ) GO  TO  110 
I GROUP  = I 


no  • oo  oooo  no  oooo  oo 
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60  TO  (1000*200*200*200*200*200*200*8000*9000*300*9200) * IGROUP 
110  CONTINUE 

THE  GROUP  CODE  IS  INCORRECT. 

120  FMT(l)  = (+6HGR0UP  ) 

130  PRINT  132*  FMT(l) 

132  FORMATt  10X  14HTHE  FOLLOWING  » A6*  18HC0DE  IS  INCORRECT.  ) 

140  PRINT  142*  ( TEST(J) *d=l,9  ) 

142  FORMAT ( 24 X 6A6*  3F12.5  ) 

150  print  152 

152  format ( iox  28hthe  program  cannot  continue.  ) 

STOP  6 


check  the  coord  code. 


200  DO  210  1=1*8 

IF C TEST(3).GT,C00RD(I)  .OR.  TEST(3) .LT.COORD(I)  ) GO  TO  210 
1COORO  = I 

GO  TO  (300*300*300*500*300.300.600,500)*  IC00RD 
210  CONTINUE 

THE  COORD  CODE  IS  INCORRECT. 

FMT(l)  = (+6HC00RD  ) 

GO  TO  130 


CHECK  THE  UNITS  CODE. 

300  DO  310  1=1*30 

IF C TEST (5). GT, UNITS (2*1)  .OR.  T£ST(5) .LT.UNITS(2* I)  ) GO  TO  310 
I UNITS  = I 

IF ( IGROUP. EQ. 10  ) GO  TO  320 
IF ( IC00RD.GE.6  ) GO  TO  400 
IF ( TEST (9)  ) 305,315*305 

305  NCARDS  = <TEST(8)-TEST(7) )/TEST(9)  ♦ 1.5 
IF C NCARDS* (1001-NCARDS)  ) 2100,2100*320 

310  CONTINUE 

the  units  code  is  incorrect. 

FMT(l)  = (+6HUNITS  ) 

GO  TO  130 

PRINT  THE  DATA, 

315  TEST(8)  = TEST (7) 

ncaros  = 1 

320  IF ( LINES. LT.LINMAX-5  ) GO  TO  330 
PRINT  322 
322  FORMAT ( 1H1  ) 

LINES  = 1 
GO  TO  340 

330  IF ( IHEAO.EQ.1  ) GO  TO  350 
PRINT  332 
332  FORMAT ( //  ) 

LINES  = LINES  * 3 

340  PRINT  342  i'* 

342  FORMAT ( 25X  73HGROUP  COORDINATE  UNITS  INITIAL 

1 FINAL  INCREMENT*  / ) 23 


oooo  oo  oooo  oo  oo 
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LINES  s LINES  ♦ 2 
IHEAD  = 1 

SET  THE  VARIABLE  FORMAT. 

350  FMT(l)  z (+6H!  24X  ) 

FMT (2)  z U6H6A6,  3) 

IF ( UNITS  d » IUnITS) , 6T.2.0E-3  ) 60  TO  360 

FMT (3)  z (+6HF12.2)) 

60  TO  370 

360  FMT (3)  z (♦6HF12.5)) 

370  PRINT  FMT , (TEST( J) , Jzl,9) 

LINES  z LINES  ♦ 1 

TEST(7)  z TEST (7)  * UNITS! 1, IUNITS) 

TcST(8)  z TEST (8)  e UNITS(1, IUNlTS) 

TEST(9)  z TEST (9)  * UNITS(1. IUNITS) 

60  TO  ( 380 • 2000 > 3000. 4000 > 5000 .60 00 » 38o»380. 360. 9100. 100) * I6R0UP 

the  codes  are  inconsistent. 

380  PRINT  382 

382  FORMAT!  10X  33HTHE  ABOVE  COOES  ARE  INCONSISTENT.  ) 

60  TO  150 


CHECK  THE  AXES  CARD. 


400  IF!  NTAPE1.EQ.0  > NTAPElzi 
IF!  IUNITS, 6T. 4 ) 60  TO  420 

IPAPER  z 600  - 500  * i ( IUNITS-D/2) 

DO  410  Jzl.3 

IF!  TEST ( J+6)  * !100.0*UNITS(1.3)-TEST(J+6)*UNITS(1, IUNITS) ) ) 
1 420.410,410 

410  CONTINUE 

60  TO  (420,420,420,4100,5100,6100.7100,420,420),  I6R0UP 

THE  AXES  CARD  IS  INCORRECT. 

420  PRINT  422 

422  FORMAT!  10X  37HTHE  FOLLOW I NO  AXES  CARD  IS  INCORRECT.  ) 

60  TO  140 
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60  TO  100 

540  IF ( IC00RD.EQ.4  ) GO  TO  550 
TEST (12)  = SPHASE(FREQ) 

GO  TO  100 

550  TEST(12)  s SLQSS(FREQ) 

GO  TO  100 


check  the  number  of  caros  in  the  profile. 


600  NCAROS  = TEST (7) 

IF ( (NCAROS. GE. 2). AND. (200. GE. NCAROS)  ) GO  TO  610 
PRINT  602.  (TEST(J) ,J=1.4) , NCaRDS 
602  FORMAT ( 10X  27HTHE  NUMBER  OF  CAROS  IN  THE  » 2(A6,A4)» 

1 110.  44h»  EXCEEDS  200  CAROS  OR  IS  LESS  THAN  2 CAROS.  ) 
GO  TO  150 

PRINT  THE  HEADING. 

610  I HEAD  = 0 

LINES  = NCAROS/2  ♦ 10  ♦ LINES 
IF ( LINES. LT.LINMAX-1  ) GO  TO  620 
PRINT  322 

LINES  = NCARDS/2  ♦ 6 
GO  TO  630 
620  PRINT  332 

630  PRINT  632.  (TEST(J) »J=1,4) 

632  FORMAT!  51X  2(A6,A4).  /.  5lX  -»  / ) 


CHECK  THE  UNITS  CODE. 

700  READ  101,  TEST(l) ,TEST(10) ,TEST(2) .TEST(ll) »TEST(3) ,TEST(12) 

J = 1 

710  DO  730  1=1,32 

IF ( TEST (J).GT.UNITS(2,I)  .OR.  TEST(J) .LT .UNITS(2»I)  ) GO  TO  730 
IF ( J.GT.l  ) GO  JO  720 
I UNITS  = I 
J = 2 
GO  TO  710 

720  IF  (J  ,GT.  2)  GO  TO  725 
JUNITS  = I 
J = 3 

IF  (IGROUP  .EQ,  7)  GO  TO  710 
GO  TO  740 
725  KUNITS  = I 
GO  TO  740 
730  CONTINUE 

THE  UNITS  CODE  IS  INCORRECT. 

FMT(l)  = (+6HUNITS  ) 

PRINT  132.  FMT(l) 

PRINT  142,  TEST ( 1 ) , TEST ( 3 ) » TEST ( 2 ) » TEST ( 4 ) 

GO  TO  150 

SET  THE  VARIABLE  FORMAT. 

740  FMT(l)  = (*6H<  2(  ) 

FMT (2)  = (+6H  115,  ) 

IF ( UNlTS(l»IUNlTS).GT.2.0E-3  ) GO  TO  750  25 


UUUUU  UOOOO  u u uuuu 
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FMT(3)  = (+6HF14.2,) 

60  TO  760 

750  FMT<3)  = (+6HF14.5* ) 

760  IF  ( I GROUP  ,E0.  7)  GO  TO  763 
LUNITS  = JUNITS 
GO  TO  765 

763  LUNITS  = KUNITS 

765  IF  ( UNITSU, LUNITS)  ,GT.  2.0E-3  ) GO  TO  770 
FMT(4)  = (♦6HF16.2#) 

GO  TO  780 

770  FMT(4)  = (+6HFI6.5.) 

780  FmT  t5)  = U6H15X)  )) 

GO  TO  (380*380*380*4200*5200*6200*7200*380*380) * IGROUP 


READ  AND  WRITE  A COMMENT  CARO. 

1000  READ  1001>  (TEST(J) > J=l>12) 

1001  FORMAT < 12A6  ) 

1002  FORMAT ( 10X  12a6  ) 

IF ( LINES. LT.LINMAX  ) GO  TO  I0O5 
PRINT  322 
LINES  = 1 

1005  PRINT  1002*  (TEST(J) >U=1#12) 

LINES  = LINES  * 1 
ICMMNT  = ICMMNT  + 1 
IF { ICMMNT, GT.l  ) GO  TO  100 
DO  1010  J=l*12 

1010  CMMNT 10)  s TEST(J) 

GO  TO  100 


CONVERT  ANO  STORE  THE  TARGET  DATA. 


IF ( IC00RD.GT.2  ) GO  TO  380 
IF ( NTAPE2.EQ.0  ) NTAPE2=1 
TARGET (1004* ICOORD)  = UNITStl* IUNITS) 
TARGET(1005>ICOORD)  = UNITS(2*IUNITS> 

N = TARGET(1010>ICOORO) 

NCARDS  = N ♦ NCARDS 
IF ( NCAROS.GT.lOOO  ) GO  TO  2100 
TARGET ( 1010* ICOORD)  = NCARDS 
N = N ♦ 1 

DO  2010  i=n>ncards 

TARGETd, ICOORD)  = TEST(7)  ♦ <I-N)*TEST(9) 
GO  TO  100 


THE  NUMBER  OF  ENTRIES  HAS  BEEN  EXCEEDEO. 
PRINT  2102*  GROUPUGROUP)  > COORO ( ICOORO) 
FORMAT ( 10X  14HTHE  NUMBER  OF  > A6>  1XA5» 
1 23HS  EXCEEDS  1000  ENTRIES,  ) 

GO  TO  150 


CONVERT  ANO  STORE  THE  SONAR  DATA 


26 
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3000  IF(  ICOORO-3  ) 3100*3200*380 

3100  SONAR (1.IC00RD)  = TEST(7) 

SONAR (2* I COORD)  = TEST(8) 
S0NAR(3,IC00R0)  = TEST (9) 

SONAR (4* ICOORO)  = UNITS(1*IUNItS) 
SONAR  IS * 1CO0R0)  = UNITS12* IUN1TS) 
GO  TO  100 


STORE  THE  SONAR  ANGLE  DATA. 

IF(  ABSCTEST(7)).GT. 1.5707  .OR.  TESTIS) ,6T. 1.5707  > GO  TO  3300 
ANGLES(1004)  = UNITS (1  * IUNITS) 

angles ( 1005)  = unitsu* iunits) 

N = ANGLES(lOlO) 

NCARDS  = N + NCARDS 

IF ( NCARDS, GT. 1000  > GO  TO  2100 

angles (loio)  = ncards 

N = N ♦ 1 

00  3210  I=N* NCARDS 

ANGLES ( I ) = TEST (7)  ♦ ( I-N) *TEsT{9) 

GO  TO  100 


THE  SONAR  ANGLES  EXCEED  THEIR  BOUNDS. 

3300  PRINT  3302 

3302  FORMATU0X50HALL  SONAR  ANGLES  MUST  LIE  BETWEEN 
GO  TO  150 


-90  AND  +90  DEG.) 


CONVERT  AND  STORE  THE  SURFACE  DATA. 


IF ( ICOORD.NE.2  ) GO  TO  380 
SURFAC(l)  s TEST (7) 

SURF AC (2)  = TEST (8) 

SURFACC3)  = TEST (9) 

SURFACE)  = UNITS  d*  IUNITS) 

SURFAC (5)  = UNITS (2* IUNITS) 

GO  TO  100 

STORE  THE  SURFACE  AXES  DATA. 

SP (207*1)  = 2. 777777777E-4/UNlTS(l, IUNITS) 
SP(207*2)  = SP{207*1) 

SP(208*1)  = TEST (7) 

SP(208*2)  = TEST (8) 

GO  TO  5100 


READ  AND  PRINT  THE  SURFACE  PROFILE. 


4200  SP(204*1)  = UNITS(1# IUNITS) 

SP(205»1)  = UNITS(2* IUNITS) 

SP(204*2)  s UNlTSd* JUNITS) 

SP(205i2)  = UNITS (2»UUNITS) 

SP(210>1)  = NCARDS 
SP(210*2)  = TEST (5) 

SURFAC (4)  = 0.0 

PRINT  4202*  SP(205*1)*SP(205*2>*SP(205,1)*SP(205*2) 

4202  FORMAT ( 2(12X  16HCARD  RANGE-*  A6*  9H  DEPTH-*  A6.  11X)» 


•p.*— --.‘■rJ., 
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REAO  4212,  (SP(N,1),SP(N»2) ,N=1» NCARDS) 

4212  FORMAT ( 2F10.5  ) 

J = ( NCANOS  ♦ 1 ) / 2 
00  4220  Lsl' J 
N s L ♦ J 

IF ( N.LE. NCARDS  ) 60  TO  4220 

PRINT  FMT » L'SP(L'1)'SP(L'2) 

60  TO  4230 

4220  PRINT  FMT'  L'SP(L' 1) »SP(L'2) » N'SP(N'l) 'SPIN. 2) 

CONVERT  AND  CHECK  THE  PROFILE. 

00  4260  Nsl.NCARUS 
SP(N.l)  = SP(N.l)  * SP(204.1) 

SP ( N ' 2 ) s SP;n,2)  * SP(204,2) 


is  the  range  increasing, 

IF ( N.EQ.l  ) GO  TO  4260 

IF ( SP(N,1).GT.SP(N-1»1)  ) GO  TO  4260 

NCAROS  = N 

GO  TO  4270 

CONTINUE 

GO  TO  100 


A CARO  IS  INCORRECT  OR  OUT  OF  ORDER. 

4270  PRINT  4272,  NCARDS 

4272  FORMAT ( 10X  llHCARD  NUMBER,  14,  13hIS  INCORRECT. 
GO  TO  150 


CONVERT  ANO  STORE  THE  BOTTOM  DATA, 


IF ( IC00R0.NE.2  ) GO  TO  380 
BOTTOM(l)  = TEST<7> 

BOTTOMU)  s TEST (8) 

BOTTOM(3)  s TEST (9) 

B0TT0M<4)  s UNITSU'IUNITS) 

B0TT0M(5)  s UNITS(2'IUNITS) 

GO  TO  100 

STORE  THE  BOTTOM  AXES  OATA. 

5100  BP(207,1)  s 2.777777777E-4/UNlTS(l,IUNlTS) 

BP(207,2)  s BP (207,1) 

BP(208, 1 ) s TEST (7) 

BP(208,2)  s TEST (8) 

60  TO  100 


REAO  AND  PRINT  THE  BOTTOM  PROFILE. 


5200  BP <204, 1 ) s UNlTS(l'IUNITS) 

BP(205,1)  s UNITS(2, I UNITS) 

6P  (204,2)  = UNITSU'JUNITS) 

BP(205,2)  s UNITS(2,UUNITS) 

BP (210,1)  s NCARDS 
BP(210'2)  s TEST (5) 

B0TT0M(4)  s 0.0 

PRINT  4202,  BP (205, 1 ) ,BP< 205,2 ) , BP (205, 1 ) ,BP( 205,2) 


»»<*«««  V<: 
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READ  4212,  (BP(N,1) »BP(N»2) 'N=X#NCAR0S> 

J = ( NCAHOS  ♦ l ) / 2 
00  5220  L=1»J 

N = L ♦ J 

IF ( N.LE.NCARDS  ) 60  TO  5220 

PRINT  FMT*  L.BP(L'l) *BP(L'2) 

GO  TO  5230 

5220  PRINT  FMT » L*BP(L*  1)  »BP',L#2) » N»BP(N#1)  *BP(N,2) 

CONVERT  AND  CHECK  THE  PROFILE. 

DO  52b0  Nzl'NCAROS 
bP(N'l)  s BP(N'l)  * BP(204,1> 

BP(N*2)  : BP(N,2)  * BP (204 .2) 

IS  THE  RANGE  INCREASING. 

IF(  N.EO.l  ) GO  TO  5260 
IF C BP(N,1).GT.BP(N-1'1)  ) GO  TO  5260 
NCARDS  = N 
GO  TO  4270 
CONTINUE 
GO  TO  100 


CONVERT  AND  STORE  THE  VELOCITY  DATA. 


6000  IF C IC00RD.NE.5  ) 60  TO  380 

TOLERA(l)  = TEST (7) 

T0LERA<2)  = TEST (8) 

TOLERAO)  s TEST (9) 

T0LERA<4)  = UNITS ( 1 » IUN ITS) 
TOLERA(5)  = UNITS(2» IUNITS) 
GO  TO  100 


6100 


STORE  THE  VELOCITY  AXES  DATA. 


VP (207  r 1 ) 
VP(207.2) 
VP(208f 1} 
VP(208'2) 
60  TO  100 


2.777777777E-4/UNlTS(l, IUNITS) 
VP(207'1) 

TEST (7) 

TEST (8) 


READ  ANO  PRINT  THE  VELOCITY  PROFILE. 


6200  VP(204'1)  = UNITSll' IUNITS) 

VP(205'1)  = UNITS(2» IUNITS) 

VP(204'2)  = UN I TS ( 1 » JUN ITS) 

VP(205'2)  = UNITS(2» JUNITS) 

VP (210 ' 1 ) s NCARDS 
VP(210.2)  = TEST (5) 

PRINT  6202'  VP(205'1)'VP(205'2)'VP(205'1)'VP(205'2) 

6202  F0RMAT(2(12X  16HCARD  DEPTH-*  A6*  lOH  VELOCITY-'  A6*  10X) * /) 

REAO  4212*  (VP(N'l)'VP(N'2)'N=l'NCAR0S) 

6210  J = ( NCARDS  ♦ 1 )/2 
00  6220  Lsl'J 
N s L ♦ J 

IF < N.LE.NCARDS  ) 60  TO  6220 

PRINT  FMT»  L*VP(L*1)'VP(L'2) 


• «V 
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GO  TO  6230 

6220  PRINT  FMT»  L*VP(L»1) »VP(L»2) » N* VP(N* l) *VP(N,2) 


CONVERT  AND  CHECK  THE  PROFILE, 

DO  6260  N=l» NCAROS 

VP(N» 1)  = VP(N,1)  * VP(204»1) 

VP(N»2)  = VP(N,2>  * VP (204, 2) 

IF(  VP(N,2).GT.1.55  .AND,  1.77.GT.VP(N,2>  ) GO  TO  6240 


THE  VELOCITY  EXCEEDS  ITS  BOUNDS, 

PRINT  6232,  N 

FORMAT ( 10X  26HVEL0CITY  PROFILE  CARD  NUMBER*  14, 
1 20H  EXCEEDS  ITS  BOUNDS,  ) 

60  TO  150 


IS  THE  DEPTH  INCREASING. 
5240  IF < N.EQ.l  ) 60  TO  6260 

IF { VP(N,1).GT.VP(N-I,1) 
NCAROS  = N 
GO  TO  4270 
6260  CONTINUE 
GO  TO  100 


GO  TO  6260 


r 


CONVERT  AND  STORE  TEMPERATURE  DATA, 


STORE  THE  THERMAL  AXES  DATA, 

TP(207,1)  = 2.777777777E-4/UNITS<l,IUNlTS) 
TP(207,2)  = TP(207,1) 

TP(207,3)  S TP(207, 1) 

TP (208,1)  = TEST  (7) 

TP(208,2)  = TEST  (0) 

TP (208,3)  s TEST  19> 

GO  TO  100 

READ  AND  PRINT  THE  THERMAL  PROFILE. 


7200  VP(204»1 

VP(205, 1 
VP(204,2 
VP (205,2 
VP(210, 1 
VP(210,2 
TP(204,l 
TP(205, 1 
TP(204»2 
TP(20b*2 
TP(204,3 
TP(205,3 
TP(210, 1 
TP(210,2 
TP(210,3 


S UNITS(1, IUNITS) 
s UNITS(2, IUNITS) 

= UNITS(1,KUNITS) 

= UNITS (2»KUNITS) 

= NCAROS 
= TEST(5) 

SUNITS(1, IUNITS) 

= UNITS(2, IUNITS) 

= UNlTSd,  JUNITS) 

= UNITS (2, JUNITS) 

= UNITS(1»9) 
s UNITS (1,10) 

= NCAROS 
= TEST (8) 
s TEST  (9) 

PRINT  7202,  TP (205,1) , TP (205, 2) *TP (205,1) , TP (205,2) 

7202  FORMAT ( 2 (6X, 16HCARD  DEPTH-, A6 , 12HTEMPERATURE- , A6, 14HSALINITY 

30 
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1-/1000. 6X)t/> 

READ  7205.  (TP(N,1) ,TP(N»2) ,TP(N, 3) »N=1, NCAROS) 

7205  FORMAT(3F10.5) 

C 

C TEST  WHICH  SALINITY  INPUT  TO  USE. 

C 

SAL  = 0.0 

00  7207  I = 1 . NC ARDS 

IF  ( TP (I .3)  ) 7208.  . . 

SAL  = SAL  ♦ TP (I, 3) 

IF  (SAL  .ST,  0.0)  60  TO  7215 

7207  CONTINUE 
GO  TO  7210 

7208  PRINT  7237,  I 
PRINT  7209 

7209  FORMAT  ( 1 0 X , 36hNEG AT I VE  SALINITY  IS  NOT  VALID  DATA.) 

GO  TO  150 

7210  DO  7212  I = l.NCARDS  ' 

TP  (1 .3)  s TPU10.3) 

7212  CONTINUE 
C 

C PRINT  THE  TEMPERATURE  PROFILE. 

C 

7215  J t ( NCAROS  ♦ 1 )/2 
DO  7220  L=1’J 
N = L ♦ J 

IF  (N  .LE.  NCAROS)  GO  TO  7220 
PRINT  7225,  L.TP(L.l) »TP(L,2) »TP(L,3) 

GO  TO  722b 

7220  PRINT  7225,  L.TP(L.l) ,TP(L,2) »TP(L,3) .N.TP(N.l) ,TP(N,2> ,TP(N,3> 

722b  F0RMAT(2(bX, I3,6X,F8.2,8X,F7.2, 11X.F7.2, 10X) ) 

722b  PRINT  7227,  TP(210,2) 

7227  FORMAT (//,45X, 12HLATITUDE  = ,Fl0.5,8H  DEGREES,//) 

C 

C CONVERT  and  CHECK  TEMPERATURE  PROFILE 

C 

7230  DO  7260  N=1»NCARDS 

IF  (JUNITS  .EG.  31)  GO  TO  7235 
TP(N»2)  S (TP(N,2)  - 32.)  * TP(2o4,2) 

7235  IF  (TP(N,2)  .GT.-3.0  .AND.  35.0  .GT.  TP(N,2))  GO  TO  7250 
C 

C THE  TEMPERATURE  EXCEEDS  NEPTUNIAN  BOUNOS. 

C 

PRINT  7237,  N 

7237  FORMAT(10X,31HTEMPERATURE  PROFILE  CARD  NUMBER, 14, 19HEXCEEDS  ITS  B 
10UN0S.) 

C 

C CHECK  IF  SALINITY  IS  WITHIN  NEPTUNIAN  BOUNDS. 

C 

7250  IF  (SAL  .GT.  0.0)  GO  TO  7254 
K s N 

IF  (K  .GT.  1 ) GO  TO  7260 

IF  (TP(2lO,3)  ,GE.  0.0  .ANO.  43.0  .GT.  TP(2l0,3))  GO  TO  7260 
GO  TO  7256 

7254  IF  (TP(N,3)  ,GE.  0.0  .AND.  43.0  .GT.  TP(210,3))  GO  TO  7260 
C 

C THE  SALINITY  EXCEEOS  NEPTUNIAN  BOUNDS. 

C 31 
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7255  PRINT  7257 
PRINT  7237,  N 
GO  TO  7260 

7256  PRINT  7257 

7257  FORMAT(10X,40HSALINITY  IS  NOT  *ItHIN  NEPTUNIAN  BOUNOS.) 
7260  CONTINUE 

CALL  BT(TP» 

C 

C PRINT  THE  HEADING 

C 

I HEAD  = 0 

LINES  s NCARDS/2  ♦ 10  ♦ LINES 
IF  CLINES  .LT.  LINMAX  - 1)  60  TO  7320 

PRINT  322 

LINES  = NCARDS/2  ♦ 6 
GO  TO  7330 
7320  PRINT  332 
7330  PRINT  7335 

7335  FORMAT  ( 5IX, X7HVEL0CITY  PROFILE  »/»51X»17H-— 

C 

C PRINT  THE  VELOCITY  PROFILE. 

C 

PRINT  6202,  VP (205,1) , VP (205, 2) ,VP( 205,1) » VP (205 ,2) 

GO  TO  6210 
C 

c 
c 

C PROCESS  CONGRATS, 

C 

8000  IF ( PRINT. LT,TEST(3)  .OR.  PRINT. GT.TEST(3)  ) GO  TO  8010 
NPRINT  s 1 
GO  TO  8020 
8010  NPRINT  = 0 

8020  FLO ( 27 , 9 , PROCES ) s NPRINT  ♦ 8 * (NT APE 1 ♦ B*NTAPE2) 

C 

C IS  THERE  AT  LEAST  ONE  COMMENT  CARD. 

IF ( ICMMNT.GT.O  ) GO  TO  8100 
PRINT  1001,  (CMMNT(J),J=1,12) 

LINES  s LINES  ♦ 1 
C 

C INITIALIZE  THE  PLOTTER  TAPE. 

8100  IF ( NTAPE1.NE.1  ) GO  TO  8200 
CALL  PLOTS ( DATA, 2048,1  ) 

CALL  PLOT (0.0, 0.0, IPAPER ) 

CALL  PL0T(5, 0,0. 0,-3) 

CALL  SYMBOL (0.0,0.0,0.14  » CMMNT (1), 90. 0,72) 

CALL  PL0T(5. 0,0. 0,-3) 

NTAPEl  s 2 
C 

C SORT  THE  TARGET  COORDINATES, 

8200  DO  8260  Jsl,2 

N s TARGET(1010,J) 

IF < N.LT.2  ) GO  TO  8260 
LSI 

8210  X s L ♦ 1 

8220  IF ( TARGET (L, J)-TARGET ( I , J)  ) 8250,8230,8240 

8230  TARGET(X,J)  s TARGET (N, J) 

N s N - 1 
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6240 

8250 

8260 

C 

C 

8270 

C 

c 

8300 

8305 

8310 

6320 

8330 

8340 

8350 


8400 

C 

C 

C 

c 

c 

9000 

9002 

9030 

9090 

C 

c 

c 

c 

c 

c 


60  TO  8250 

TARGET(N*1. J)  = TARGET! I rJ) 
TARGET ( I . J)  = TARGET (L»d) 
TARGET (L,J)  = TARGET (N4l»J) 
I = X ♦ 1 

IF ( I.LE.N  ) GO  TO  8220 
L = L ♦ 1 

IF < L.LT.N  > GO  TO  8210 
TARGET ( lOlO’ J)  = N 
CONTINUE 
CALL  BTPLOT(TP) 


SET  THE  MAXIMUM  REVERSAL  INCREMENT. 

IF ( REVMAX(2).GT.REVMAX(1)  ) GO  To  8270 
REVMAX (3)  = 0.0 
GO  TO  8300 
N = TARGET(1010»1) 

REVMAX (3)  = (TARGET(N»1)-TARGET(1»1))/(REVMAX(2)-REVMAXU)J 

SORT  THE  SONAR  ANGLES. 

N s ANGLES (1010) 

IF ( N.LT.l  ) GO  TO  8 400 
DO  8305  Lsl'N 

IF ( ABS ( ANGLES (L)).GT.l. OE-4  ) GO  TO  8305 

ANGLES(L)  = 0.0 

CONTINUE 

IF ( N.LT.2  ) GO  TO  8400 
L = 1 
I = L ♦ 1 

IF C ANGLES (L) -ANGLES ( I)  ) 8350 » 8330 » 8340 

ANGLES! I)  r ANGLES (N) 

N = N - 1 
GO  TO  8350 

ANGLES (N+l)  = ANGLES! I) 

ANGLES! I)  = ANGLES (L) 

ANGLES  (L)  = ANGLES  (Nd) 

I = I ♦ 1 

IF!  I.LE.N  ) GO  TO  8320 
L = L ♦ 1 

IF!  L.LT.N  ) GO  TO  8310 
ANGLES! 1010)  s N 
RETURN  * 


TERMINATE  CONGRATS. 

PRINT  9002 

FORMAT!  1H1 • 9X  2BHC0NGRATS  HAS  BEEN  COMPLETED.  » /»  1H1  ) 
IF!  NTAPE2.EQ.0  ) GO  TO  9090 
END  FILE  2 
STOP  5 


READ  MISCELLANEOUS  DATA. 
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C READ  THE  FREQUENCY. 

9100  TARGET ( 1008*3)  = TEST(7) 

IF ( (IUNlTS-27)*(IUNITS-28)  ) 380,9110*380 

9110  TAR6ET11008.1)  = ATTEN(TARGET(l008,3) > 

PRINT  9112.  TARGET (1008.1) 

9112  FORMAT ( 24X,  12HATTENUATI0N  . l2X  6HDB/KYP.  6X  F12.S  ) 

LINES  = LINES  ♦ 1 
GO  TO  100 

READ  THE  MAXIMUM  NUMBER  OF  REVERSALS. 

9200  REVMAX(l)  = TEST(7) 

REVMAXC2)  s MAX(TEST<7) ,TEST(8) ) 

REVMAX(4)  s TEST (91 
TEST(8)  = REVMAX(2) 

I UNITS  = 15 

IF ( REVMAX(2).LT. 201,0  ) GO  TO  320 
PRINT  9202 

9202  FORMAT C 10X45HTHE  MAXIMUM  NUMBER  OF  REVERSALS  EXCEEDS  200.0  ) 
GO  TO  140 
END* 
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FUNCTION  BLOSS  AND  SLOSS  FOR  CONGRATS. 


CHIC  WEINBERG 


READ  AND  PRINT  THE  LOSS  TABLE. 


FUNCTION  BL0SS<  freq,test»por,dbmin»dbmax  ) 

DOUBLE  PRECISION  R,Z,RB,ZB»C0S28»RA02B*S1N2B 
DIMENSION  DB<92,4)/276*0. 0,92*3. 14159265/ 

DIMENSION  TYPE(10)/60H  AMOS  MGSl  MGS2  MGS3  MGS4  MGS5 

1NEL  NUWC  OTHERS 
PRINT  2 

2 FORMATt  54X  11HB0TT0M  LOSS,  /,  54X  11H / ) 

J = 1 

DO  10  1=1,10 

IF ( TEST.GT.TYPE(I)  .OR,  TEST.lT.TYPE( I)  ) GO  TO  10 
ITYPE  =1-1 

IF!  -ITYPE  .LE.  0 ) GO  TO  100 

CALL  BTMJ.0S!  FREQ, DB> ITYPE, POR,DBMiN,DbMAX,*110  ) 

10  CONTINUE 

ENTRY  SLOSS ( FREQ  } 

PRINT  12 

12  FORMATC  54X  12HSURFACE  LOSS,  /,  54X  12h— — — — , / ) 

J = 2 
GO  TO  100 

ENTRY  BPHASE(FREQ) 

PRINT  22 

22  FORMAT!  54X  12HB0TT0M  PHASE,  /,  54x  12H— — — — — , / ) 

J = 3 
GO  TO  100 

ENTRY  SPHASE(FREQ) 

PRINT  32 

32  FORMAT!  54X  13HSURFAC  PHASE,  /»  54X  13H- — — , / ) 

J = 4 

100  READ  101,  !DB!I,U) ,1=1,91) 

101  FORMAT!  10F5.1  ) 

110  D8!92, J)  = 08(91, J) 

IF!  U.GT.2  ) GO  TO  150 
PRINT  112 


112 

FORMAT!  5I24H 

ANGLE 

LOSS 

), 

1 

/,  5I24H 

DEG 

db 

), 

/ 

) 

GO  TO  160 

150 

PRINT  152 

152 

FORMAT!  5124H 

ANGLE 

phase 

), 

1 

/,  5124H 

OEG 

DEG 

), 

/ 

) 

160  DO  170  1=1,18 

11  = I - 1 

12  = II  ♦ 19 

13  s 12  ♦ 18 

14  = 13  ♦ 18 

15  = 14  ♦ 18 

170  PRINT  172,  U,DB!U+1,U),  I2,DB!I2*1»J) , I3,DB(I3+1,J) , 
1 14,DB!I4+1, J) , I5,DB( I5+1, J) 
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172  FORMAT!  I8,F11.2,  I13.FU.2,  U3.F11.2,  U3»F11.2»  113, FU. 2 ) 

II  = 16 

PRINT  172,  Il»0B(U+l,J) 

IF  ( I TYPE  .GT,  0 ) PRINT  174,  TYPE(lTYPE+l) 

174  FORMAT(//»40X, ‘METHOD  OF  BOTTOM  LOSS  COMPUTATIONS  *,A6) 

IF  ( POR  ,GT.  0.0  ) PRINT  175,  pOR 
17t>  FORMAT  ( 52X,  ‘POROSITY  = »,F4.2) 

IF  (OBMIn  .GT.  0,0  .OR,  DBMIN  ,LT.  0.0  ) PRINT  176,  OBMIN 

176  FORMAT  ( 41X, ‘MINIMUM  LIMIT  FOR  BOTTOM  LOSS  = «,F6.2) 

IF  ( DBMAX  ,GT.  0.0  .AND.  DBMAX  .LT.  500.0  ) PRINT  177,  OBMAX 

177  FORMAT  ( 41X, 'MAXIMUM  LIMIT  FOR  BOTTOM  LOSS  = *»F6.2> 

FREQUE  = FREQ 

IF ( J.LE.2  ) QO  TO  190 

CONVERT  THE  PHASE  TO  RADIANS. 

DO  160  1=1,92 

ISO  DB(I.J)  s 08( I, J)  * 1.74532925E-2 

190  RETURN  r 


USE  THE  LOSS  TABLE. 


ENTRY  BOTLOS(R,Z,RB, ZB, PHASE) 

J = i 
GO  TO  200 

ENTRY  SURLOS(R,Z,RB, ZB, PHASE) 

J = 2 

200  IF ( R+RB+Z+ZB  ) 220,210,220 

210  ANGLE  = 90.0 
GO  TO  230 

220  ANGLE  = 57.295780+0  * ATAN( (R*ZB-Z*RB> / (R+RB+Z+ZB) ) 

230  I = ABS(ANGLE)  + 1.0 
ANGLEI  =1-1 

BLOSS  = OB(I.J)  + (DB(I+1,J)-DB(I,J) > * ( ABS (ANGLE) -ANGLEI ) 

PHASE  = DB(I,J+2)  + (DB ( 1+1 , J+2) -DB ( I , J+2) ) * (ABS (ANGLE) -ANGLE I) 


compute  the  direction  numbers  op  the  reflected  ray. 


RAD2B  = RB«RB  + ZB*ZB 

C0S2B  = (RB+ZB)«(RB-ZB)  / RAD2B 

S1N2B  = 2. 0D+0*RB*ZB/RAD2B 

RAD2B  = C0S2B+R  + SIN2B+Z 

Z = SIN2B+R  - C0S2B+Z 

R = RA02B 

RETURN 

ENO 
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SUBROUTINE  BTMLOS  ( FREQ  > DB » I TYPE , PoR > DBMIN > DBMAX > $ ) 

c 

DIMENSION  08(91) 

DIMENSION  FR(b)/0.1>0.3>0.5>1.0>1.5>2.5>4,0>6.5/ 

DIMENSION  FVAL(8)/0.157>0.45>0.67>1,0>1,18>1.43>1.7>2.0/ 
GO  TO  (100>1000>2000>3000>4000.5000>6000>7000>9000)> I TYPE 


BOTTOM  LOSS  FROM  AMOS  DATA 


100  FREQUE  = FREQ 

FKQLG2  S LOG (FREQUE ) /LOG ( 2 , 0 ) 

GRAZING  ANGLES  ARE  0 THROUGH  14. 

F = FRQLG2  '' 

FSQ  = F*F 

F CUBE  = FSQ*F 

F4P4R  = FCUBE*F 

Ob(l)  = 0,0 

DO  110  I = 1 > 14 

THETA  = I 

Th  = 0,2*(THETA  - 10.0) 

110  08(141)  = 2,9o026  ♦ 1.26990*F  + 1.60430*TH  4 0.25257*FSQ 

1 ♦ 0,66161*F*TH  4 0,06072*TH*TH  + 0 . 02000*FCUBE 

2 4 0.09833*FSQ*TH  - 0 . 01888*F*TH*TH  - 0 . 0Z619*TH**3 

3 - 0.00303*F4PWR  4-  0 ,00778*FCuBE*TH  - 0.00833*FSQ*TH*TH 

4 - 0.01250*F*TH**3  - 0.01310*TH**4 

TEST  FOR  FREQUENCY  GREATER  THAN  1.0  KC. 

IF  (FREQUE  ,GT,  1.0)  GO  TO  200 
FREQUENCY  IS  NOT  GREATER  THAN  1.0  KC 


F = 2.0*pRQLG2  4 3,0 
FSQ  = F*F 
FCUBE  s FSQeF 
F4PWR  = FCUBE*F 

GRAZING  ANGLES  ARE  15  THROUGH  39 


~f 

■ I 


DO  120  I=15>39 
THETA  = I 

Th  s 0,2* (THETA  - 30.0) 

120  06(141)  = 4.35024  * 1.1709l*F  * 0.46674*TH  4 0.0S272*FSQ 

1 4 0,l7115*F*TH  - 0.07547*TH*TH  - 0.003b6*FCUBE 

2 4 0.01295*FSQ*TH  - 0 ,0167B*F*TH*TH  - 0.00555*TH**3 

3 - 0.00003*F4PWR  - 0.00193*FCUBE*TH  - 0.00074*FSG*TH-»TH 

4 - 0.00224*F*TH**3 


GRAZING  ANGLES  ARE  40  THROUGH  60 
DO  130  1=40*60 

THETA  si  37 
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TH  = Ot 2* (THETA  - 50.0) 

150  0b( 1+1)  = 4.64320  ♦ 1.46858*F  - 0.11**57*TH  ♦ 0.11036*FSQ 

1 - 0.0l238*F*TH  - 0,00953*TH*tH  - 0.01000*FCUBE 

2 - 0.00375*FSQ*TH  - 0.00929*F*TH*TH  ♦ 0.00833*TH**3 

3 - 0. 00042*FCUBE*TH  - 0.00268*FSQ*TH*TH  ♦ .00830*F*TH**3 

4 -0.00208*TH**4 

GRAZING  ANGLES  ARE  61  THROUGH  90 

F = FRQLG2 
DO  140  1=61*90 

140  Dd(X+l)  = OB ( I ) - 0.085667  - 0.025434*F 
GO  TO  9000 

FREQUENCY  IS  GREATER  THAN  1.0  KC. 


200  F = FRQLG2  - 1 
FSQ  = F*F 
FCUBE  = FSQ*F 
F4PWR  = FCUBE*F 

GRAZING  ANGLES  ARE  15  THROUGH  39 

00  220  1=15*39 
THETA  = 1 

TH  = 0.2*(TH£TA  - 30.0) 

220  DB(I+1)  = 12.25083  ♦ 3.88808*F  ♦ 1.148o2*TH  ♦ 0.16312*FSQ 

1 ♦ 0 . 04386*F*TH  - 0.1&835*Th* TH  - 0.02857*FCUBE 

2 - 0.08214*FSQ*TH  - 0 . 05952*F*TH*TH  ♦ 0.00667*TH**3 

3 ♦ 0.0l071*F4PWR  - 0.00l79*FCuBE*TH  ♦ 0.00476*FSQ*TH*TH 

4 ♦ 0.01389*F*TH**3  - Q.00348*TH**4 

GRAZING  ANGLES  ARE  40  THROUGH  60 

DO  230  1=40*60 
THETA  = I 

TH  = 0.2* (THETA  - 50.0) 

230  DB(I+1)  = 13.59067  ♦ 3.74190*F  - 0.25928*TH  ♦ 0.13266*FSQ 

1 - 0.03764*F*TH  - 0.068l6*TH*TH  - 0.01333*FCUBE 

2 ♦ 0.01214*FSQ*TH  - 0.00959*F*TH*TH  ♦ 0.005C0*TH**3 

3 - 0.02000*F4PWR  - 0.00167*FCUBE*TH  ♦ 0.00867*FSQ*TH*TH 

4 - 0.00167*F*TH**3  ♦ 0.00083*TH**4 

GRAZING  ANGLES  ARE  61  THROUGH  90 


F = FRQLG2 
00  240  1=61*90 

240  0B( 1+1)  =08(1)  - 0.126002  - 0.000333*F  ♦ 0.034002/F 

GO  TO  90Q0 


BOTTOM  LOSS  FOR  ACOUSTIC  PROVINCES.  MGS  1-5. 


BOTTOM  LOSS  FOR  ACOUSTIC  PROVINCE  MSS  1. 
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00  1020  1=0*90 
Th  = I 

06(1+1)  = 2.2435*LOG(0.1260*TH  ♦ 1.496) 
GO  TO  9000 


BOTTOM  LOSS  FOR  ACOUSTIC  PROVINCE  MGS  2. 


00  2020  1=0*90  . 

TH  = I 

DUU+l)  = 3.4315*LOG(0,1056*TH  4 2.842) 
GO  TO  9000 


bottom  loss  for  acoustic  PROVINCE  MGS  3. 


DO  3020  1=0*90 
TH  = I 

06(1+1)  = 2.49l0*LOG(0.S864*TH  ♦ 10.526) 
GO  TO  9000 


BOTTOM  LOSS  FOR  ACOUSTIC  PROVINCE  MGS  4. 


00  4020  1=0*90 
TH  = I 

06(1+1)  = 2.8377*L0G(1,8754*TH  ♦ 15.68s) 
GO  TO  9000 


BOTTOM  LOSS  FOR  ACOUSTIC  PROVINCE  MGS  5. 


00  5020  1=0*90 
TH  = I 

Dd(IU)  = 2.4036*LOG(20.5760*TH  4 82.440) 
GO  TO  9000 


NAVAL  ELECTRONIC  LABORATORY  EQUATIONS. 
LOSS  DEPENDENT  ON  POROSITY  OF  BOTTOM. 


FREQUE  = FREQ 
F = 10, 0*LOG10 (FREQUE) 

POROS  = POR 
DO  6020  1=0*90 
TH  = I 

08(1+1)  = (6.0  ♦ 22.0*(POROS  - 0.27)  ♦ F)  * (TANH( (POROS/0,24)  • 
( TH/57 . 295779 ) ) ♦ (1.0  - POROS/0.24J/12.5  * 
TH*TH/8100,0> 

GO  TO  9000 
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naval  unoersea  warfare  center  equations, 
loss  dependent  on  porosity  of  bottom. 


FREQUE  = FREQ 

IF  ( FREQUE  - 0.1  ) 7030.  7030. 

IF  ( FREQUE  - 6.5  ) . , 7040 

DO  7020  1=1.8 

IF  ( (FREQUE  - FR( I ) ) *(FR( 1+1)  - FREQUE)  ) 7020.  . 

FUNU  = FVAL(I)  ♦ (FREQUE  - FR(I)  )*(FVAL(H-1>  - FVAL(I) )/(FR(I+l) 
- FR(I)  ) 


FUNU  = FVAL(I)  ♦ (FREQUE  - FR(I) )*(FVAL(I*1> 
- FR(I)  ) 

60  TO  7100 

CONTINUE  > 

FUNU  = FVAL(l) 

60  TO  7100 
FUNU  = FVAL(8) 

60  TO  7100 
POROS  = POR 

FRQFUN  = (3.7  ♦ 17.5MP0R0S  - 0.27)  )*FUNU 
PORTAN  = P0R0S/13, 75104 
SUPER  = 1 ,5/POROS 

PORDIV  = (1  - POROS/ 0. 27 )/10 1250.0 
DO  7120  1=0.90 
TH  = I 

Do(I+l)  = FRQFUN*(TANH(  (PORTAN*TH)**SuPER  ) 
60  TO  9000 


POROI V*TH*TH  ) 


CHECK  LOSS  A6AINST  MINIMUM  AND  MAXIMUM 


9000  DO  9010  1=1.91 

IF  ( D8 ( I ) ,LT.  DBMIN  ) DB(I)  = DBMIN 
IF  ( DB ( I ) ,6T.  DBMAX  ) DB(I)  = DBMAX 
9010  CONTINUE 
9900  RETURN  7 
ENO 


LO 


